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Description 

[0001] The present invention relates to image processing, for example to contrast enhancement of images by his- 
togram manipulation. 

5 [0002] A digital image is an array, usually a rectangular matrix, of pixels. Each pixel is one picture element and is a 
digital quantity that is a value that represents some property of the image at the location in the array corresponding to 
a particular location in the image. Typically, in continuous tone black and white images the pixel values represent a 
grey scale value. 

[0003] Pixel values for an image have to conform to a specified range. For example, each array element may be one 
10 byte, i.e., eight bits. In that example, the pixel values must range from 0 to 255. In a grey scale image a 255 may 
represent absolute white and 0 total black. 

[0004] Colour images consist of three colour planes, generally corresponding to red, green, and blue (RGB). For a 
particular pixel, there is one value for each of these colour planes, i.e., a value representing the red component, a 
value representing the green component, and a value representing the blue component. By varying the intensity of 

15 these three components, all colours in the colour spectrum may be created. 

[0005] Output devices, such as printers and displays, also have particular ranges for pixel values. A particular device 
may be configured to accept and output eight-bit image data, i.e., pixel values in the range 0 to 255. 
[0006] However, many images do not have pixel values that make effective use of the full dynamic range of pixel 
values available on an output device. For example, in the eight-bit case, a particular image may in its digital form only 

20 contain pixel values ranging from 100 to 150, i.e., the pixels fall somewhere in the middle of the grey scale. Similarly, 
an eight-bit colour image may also have RGB values that fall within a range somewhere in the middle of the range 
available for the output device. The result in either case is that the output is relatively dull in appearance. 
[0007] The visual appearance of an image can often be improved by remapping the pixel values to take advantage 
of the full range of possible outputs. That procedure is called contrast enhancement. 

25 [0008] There are several prior art contrast enhancement techniques. These techniques are especially common for 
monochrome (e.g., grey scale) images, such as from medical (e.g., cat-scans, x-rays, and ultrasound) and radar sourc- 
es. In such applications, contrast enhancement can be used to find image details that would otherwise be difficult to 
discern. 

[0009] Contrast enhancement techniques are often based on histogram equalisation. In histogram equalisation a 

30 histogram of grey level distribution of the image is constructed. Figure 1 is an example of such a histogram. A histogram 
is a one dimensional array with an array element corresponding to each value in the range of pixel values. Each 
histogram element contains a count of the number of pixels that has the particular pixel value corresponding to that 
element. In histogram equalisation, the pixel values in the image are altered to make the distribution of grey level values 
as uniform as possible. Histogram equalisation is described in A.K. Jain, Fundamentals of Digital Image Processing, 

35 Prentice-Hall, Inc., 1989, pp. 241-243. 

[001 0] A variation of histogram equalisation, known as adaptive histogram equalisation or local area histogram equal- 
isation, uses a sliding window to define an image region for each pixel. The histogram of that region is then equalised 
to determine the output value for the pixel. Adaptive histogram equalisation is described in John B. Zimmerman, et al., 
An evaluation of the effectiveness of adaptive histogram equalisation for contrast enhancement, IEEE Trans. On Med- 

40 jcal Imaging, 7(4):304-312, December 1988. The adaptive histogram equalisation procedure is computationally very 
expensive because a separate histogram is constructed for each image pixel. A number of region-based variants of 
adaptive histogram equalisation have been proposed in Stephen Pizer et al., Adaptive histogram equalisation and its 
variations, Computer Vision, Graphics, and Image Processing, 39(3):355-368, September 1987; S.S.Y. Lau, Global 
image enhancement using local information, Electronics Letters, 30(2) :1 22-1 23, January 1994; J. A. stark and W.J. 

45 Fitzgerald, Model-based adaptive histogram equalisation, Signal Processing, 37:193-200, 1994. These techniques 
give comparable results to adaptive histogram equalisation while reducing computation. 

[0011] Deficiencies common to adaptive histogram variants are that although some variants reduce computation 
requirements substantially relative to full-blown adaptive histogram equalisation, the required computation is still sig- 
nificant enough to be impractical in many applications. 

50 [0012] A further deficiency of adaptive histogram variants is that because these schemes adjust image pixel values 
locally, the processed image will not retain the relative lightness of pixels from the original image. For example, some 
details in a shadow region may be lightened to the point where these details are lighter than sunlit parts of the same 
image. The uneven lightening/darkening of an image leads to unacceptable visual appearance in photographic images. 
[0013] A generalisation of histogram equalisation, known as histogram specification, adjusts the data so the output 

55 histogram has a distribution close to some desired distribution. Histogram specification is described in Jain, supra, 
page 243. 

[0014] Colour image enhancement is considerably more complicated. Colour images consist of three colour planes, 
generally corresponding to red, green, and blue (RGB). One simple colour image enhancement technique (described 



2 



EP0 772 158 B1 



in RE. Trahanias and A.N. Venetsanopoulos, Colour image enhancement through 3-d histogram equalisation, 11th 
IAPR Conference on Pattern Recognition, Conference C: Image, Speech, and Signal Analysis, pp. 545-548, September 
1992) is to perform histogram equalisation on each colour independently. However, that technique can cause large 
colour shifts and other undesirable artifacts in the image. 

[0015] One histogram-based contrast enhancement technique that attempts to avoid such artifacts is to construct a 
three dimensional histogram and perform modifications in all three dimensions jointly. That technique is described in 
Trahanias and venetsanopoulos, above, and in Phillip A. Misna and Jeffrey J. Rodriguez, Explosion of multidimensional 
image histograms, IEEE International Conference on Image Processing, Vol III, pp. 958-962, Austin, Texas, November 
1994. Unfortunately, the required computation for the technique of joint modification of a three-dimensional histogram 
is at least of the order of the number of occupied histogram bins, which can number in excess of a hundred thousand 
for a moderately large image. This approach is therefore computationally intensive and hence not very practical. 
[0016] Other approaches, including the preferred embodiment of the present invention, transform the image data 
into some luminance-chrominance colour space and perform the enhancement in this domain. Some techniques adjust 
only the luminance histogram, while others adjust both luminance and chrominance histograms. Straightforward his- 
togram equalisation on the luminance component can cause undesirable artifacts in the output image. Specifically, 
large areas of approximately equal luminance can suffer from contouring artifacts, while an especially dark (light) 
background can cause the foreground objects to become too light (dark). Some of the more complicated enhancement 
procedures require that the image data undergo a non-linear transformation before processing, which increases the 
computational complexity. Histogram equalisation of luminance or chrominance histograms is described in Robin N. 
Strickland, et al., Digital colour image enhancement based on the saturation component, Optical Engineering, 26(7): 
609-61 6, July 1 987, and in Ilia M. Bockstein, Colour equalisation method and its application to colour image processing, 
Journal of the Optical Society of America A, 3(5):735-737, May 1 986. 

[0017] Histogram manipulation has also been used for purposes other than contrast enhancement. For example, for 
monochrome images, multi-thresholding techniques have been developed for dividing a histogram into clusters. The 
resulting clusters are then used to segment the image into regions corresponding to various objects and backgrounds. 
Such clustering techniques are described in P.K. Sahoo, et al., A survey of thresholding techniques, Computer Vision, 
Graphics, and Image Processing, 41 (2):233-260, February 1 988. If the histogram contains flat regions of low amplitude 
between higher peaks, for example, a threshold value anywhere in the flat region would give nearly the same segmen- 
tation of the image. For the purpose of contrast enhancement, it is desirable to differentiate between regions of high 
activity and regions with little activity. Consequently, it is desirable to divide the histogram such that flat regions of low 
amplitude form separate clusters. Therefore, the multi-thresholding clustering techniques are not suited to colour con- 
trast enhancement. 

[0018] The present invention seeks to provide improved image processing. 

[0019] According to an aspect of the present invention, there is provided a method of operating a computer to produce 
contrast enhanced images as specified in claim 1. 

[0020] According to another aspect of the present invention, there is provided an image processing system as spec- 
ified in claim 8. 

[0021] It is thus possible to provide contrast enhancement which is not computationally intensive and which does 
not introduce undesirable artifacts in the output image. 

[0022] A first embodiment provides a method of operating a computer to produce contrast enhanced digital images. 
The method commences with the step of producing a histogram having a first axis corresponding to a measurable 
property (e.g., luminance) and a second axis corresponding to a count of pixels having a particular value for the meas- 
urable property. This histogram is divided into clusters and histogram equalisation or stretching is performed on each 
cluster thereby producing a modified histogram. The modified histogram is used to adjust the value of said first meas- 
urable property in said digital form, thereby producing a contrast enhanced image. 

[0023] In one embodiment, the histogram is divided into clusters using a pattern matching technique. For example, 
patterns in the histogram that resemble gaussian distributions and patterns that resemble uniform distributions are 
separated into individual clusters. 

[0024] In another embodiment, subsequent to dividing the histogram into clusters, the cluster boundaries are rem- 
apped proportionally to the cluster's range, the number of pixels in the range, the total number of pixels in the image, 
and the histogram extent. 

[0025] In another embodiment, some other property (e.g., chrominance) of the image is used to limit the values of 
the modified histogram. 

[0026] Another embodiment provides a computer system for image processing having a histogram producing unit 
operable to produce a histogram data structure in a memory on the computer system. The computer system further 
contains a clusterizer connected to the histogram data structure and is operable to partition the histogram data structure 
into clusters. A histogram cluster boundary adjuster is connected to the histogram data structure and is operable to 
adjust boundaries of each cluster in said data structure such that said boundaries conform to a function having the 
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parameters cluster width (RO), a count of pixels in the cluster (Ni), a count of pixels in the image (Nt), the original 
histogram extent (EO), and a desired histogram extent (Ed). A histogram adjuster is also connected to said histogram 
data structure and is operable to modify the histogram for each such cluster according to a desired characteristic. 
[0027] A preferred embodiment provides apparatus for improving the visual appearance of digitised images stored 

5 in an image data structure by contrast enhancement, comprising a histogram data structure; a histogram generator 
connected to the histogram data structure and operable to store in the histogram data structure histogram data for a 
property of the images; a histogram clusterizer connected to the histogram data structure and operable to define clusters 
in the histogram, in which each such cluster has a range; a histogram remapper connected to the histogram data 
structure and operable to adjust the range of each such cluster; and a histogram adjuster connected to the histogram 

10 data structure and operable to adjust the histogram of each cluster. 

[0028] Preferably, the histogram adjuster is a histogram equaliser or stretcher. 

[0029] The apparatus may include an image transformer connected to the image data structure and operable to 
transform the image from a first data format to a second data format having the first property. 

[0030] The preferred converted image format is selected from the set of Y-Cb-Cr, CIE L*a*b*, CIE XYZ, and Y-U-V. 
15 [0031] An embodiment of the present invention is described below, by way of example only, with reference to the 
accompanying drawings, in which: 

Figure 1 is an exemplary histogram showing the number of pixels of an image having each valid grey level value. 

20 Figure 2 is a block diagram of an embodiment of computer system with enhanced contrast of digitised images. 

Figure 3 is a block diagram of the process flow of an image processing system. 

Figure 4 is a block diagram of an embodiment of image processing system having a contrast enhancer. 

25 

Figure 5 shows a block diagram of the operation of a contrast enhancer of the system of Figure 4. 
Figure 6 is an exemplary luminance histogram showing two peaks separated by a flat-region. 
30 Figure 7 graphically shows the division of a histogram into several clusters. 

Figure 8 is a flow chart of the preferred clusterizing procedure. 
Figures 9 a is an exemplary unmodified luminance histogram. 

35 

Figure 9b, 9c 1 and 9d are luminance histograms modified by the preferred contrast enhancer using different pa- 
rameter values. 

[0032] The preferred embodiment may be practised on a general purpose computer system, on a special purpose 
40 image processing system, in software used for image processing, in device drivers and other applications involving 
digital image manipulation. For illustrative purposes, the preferred embodiment is described below in the context of a 
device driver loaded into a general purpose computer for execution when an image is output on a particular device 
controlled by the device driver. 

45 A. Computer system description 

[0033] Figure 2 is a block diagram of a computer system 100 in which the preferred processing system is used to 
enhance the contrast of digitised images. The computer system 1 00 consists of a computer 1 01 connected to several 
input-output devices, e.g., a scanner 1 1 3, a printer 1 1 5, a monitor 1 1 7, and one or more disk drives 1 1 9. The computer 
50 101 may be connected to other input devices 121. These other input devices 121 may include one or more of the 
following: frame grabbers, CCD cameras, CD ROM drives, digital cameras, medical imaging devices such as CAT- 
scan, X-ray, and ultrasound. 

[0034] The computer 1 01 consists of a central processing unit (CPU) 1 03 and a random access memory (RAM) 1 05. 
The computer 1 01 may further contain a read only memory (ROM) 1 07. These devices are connected to one another 
55 via bus 1 09. The computer 1 01 may also contain input-output control units 111a through 1 1 1 f connected to the input- 
output devices 113 through 119. 

[0035] In the operation as a digital image processing system, a digitised image is input to the computer 101. The 
source of the image may be a scanner 1 1 3. The source of the digitised image may also be another source of digitised 



4 



EP0 772 158 B1 



images, such as, photo cd, a digital camera (e.g., a CCD (charged coupled devices) camera), or a frame grabber. 
[0036] Furthermore, the computer system 1 00 may be connected via a network 1 23 to other computer systems 1 00' 
and 100", etc. The digitised images may then be down-loaded onto system 100 from one or more of these other 
computer systems 100' and 100". 

5 [0037] An image may be scanned directly into the RAM 1 05 of the computer 1 01 or an image 1 09 may be stored in 
secondary storage, such as, the disk drive 1 1 9. Depending on its size, the disk drive 1 23 can store a large number of 
digitised images. The disk drive 123 also contains one or more output drivers 125. There is an output driver 125 
associated with each output device 1 1 5. The output devices are may be printers, e.g., colour ink jet printers or colour 
laser printers, such as those manufactured by the Hewlett-Packard Company, Palo Alto, California. The output devices 

10 may also include dye sublimation printers, e.g., the Kodak Colorease printer. 

[0038] When the operating system software (not shown) of the computer 1 01 encounters a command to output an 
image 1 09 on one of the output devices 115, it loads the appropriate output driver 1 25 into the RAM 1 05. 
[0039] In an alternative embodiment, the output driver is stored in ROM 1 07. In such embodiments, the step of loading 
the output driver is not necessary. 

15 [0040] Figure 3 is a block diagram of the process flow of an image processing system. An image of an object 301 is 
captured by an imaging device 303. A digitiser 305 converts the image into a digital representation. The digital version 
of the image is then stored 307 into a storage device, e.g., the RAM 1 05. Subsequent to storing, the image is retrieved 
from the storage device and processed 309. Processing 309 may either be done by the CPU 103, or by processors 
located in the input devices or output devices. Processing 309 may further be partitioned between input, output, and 

20 main computer system processors. Finally, the processed image is output 311 onto some form of output device, e.g., 
printer 115. 

[0041] It should be noted that the preferred contrast enhancer may be employed at any stage in the process flow of 
Figure 3. 

25 B. Device Drivers 

[0042] Figure 4 is a block diagram of an image processing system having a contrast enhancer according to the 
preferred embodiment. Figure 4 shows a device driver 125 loaded into the RAM 105. The output driver 125 consists 
of a contrast enhancer 205 and an output control 203. The output control 203 accepts as input digitised images and 

30 outputs control instructions to the output device 115. For a printer driver, the printer driver (output driver) includes 
instructions to control the amount and placement of ink of various colours on the output medium. 
[0043] The contrast enhancer 205 is connected to a data structure for storing a histogram 211 of some property of 
the image 109. In one embodiment, the histogram property is luminance, i.e., the histogram 211 contains an entry for 
each possible luminance value for the particular output device 1 1 5 and each such luminance value entry is the number 

35 of pixels that has that particular luminance value. 

[0044] The contrast enhancer 205 contains several modules. A first module is an image transformer 207. The image 
transformer 207 accepts as input an image 109 in one format and transforms it into an image having at least one 
property used for contrast enhancement. In one embodiment, luminance is used for contrast enhancement, thus, the 
image transformer 207 transforms the input image 109 into a representation having a luminance component, e.g., 

40 YCrCb or CIELab colour co-ordinate systems. The functionality of the image transformer 207 is described in greater 
detail below in conjunction with Figure 5. 

[0045] The contrast enhancer 205 contains a histogram generator 209. The histogram generator 209 counts the 
number of pixels having each value in the range of values allowable for the particular property used for contrast ad- 
justment of the image. These pixel counts are stored in a one dimensional histogram array 211 in the RAM 1 25. In one 
45 embodiment, luminance is used for contrast enhancement. Thus, in that embodiment, the histogram contains counts 
of pixels each having a valid luminance value. 

[0046] The contrast enhancer 205 further contains a clusterizer 213 for determining clusters in the histogram 211 . 
The operation of the clusterizer 213 is described in greater detail below in the section entitled "Clusterizer". 
[0047] Having established clusters in the histogram 21 1 , the contrast enhancer 205 calls upon a cluster adjuster 215 
50 to adjust the boundaries of the various clusters. 

[0048] Histogram equalisation or histogram stretching is then performed on the adjusted clusters by a histogram 
equaliser/stretcher 219. And using the output histogram, the pixel values are remapped by a luminance histogram 
remapper 217. 

[0049] A data checker 221 verifies that output values remain within allowed bounds for the particular colour space 
55 used by this embodiment. 

[0050] The various modules of the contrast enhancer 205 are controlled by a control program 223. The control pro- 
gram 223 implements the flow-diagram of Figure 5. 
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C. Contrast enhancer 

[0051] Figure 5 shows a block diagram of the operation of the preferred embodiment of contrast enhancer 205. The 
contrast enhancer 205 first divides the luminance histogram into clusters. The contrast enhancer 205 next performs 
5 histogram equalisation or stretching on each cluster separately, and the entire cluster is remapped to a new luminance 
region based on the cluster width, the number of pixels in the cluster, and the original histogram extent. 
[0052] The chrominance information is used to limit the modified luminance values prior to remapping the pixel values 
to avoid colour shifts due to saturation. 

10 D. Colour Transformation 

[0053] In a first embodiment, the contrast enhancement is based on luminance values. Therefore, in that embodi- 
ment, first, the contrast enhancer 205 converts the RGB image data of an input image 1 09 to a luminance-chrominance 
representation, step 501. 

15 [0054] In general terms, the preferred embodiment is equally applicable to all colour image representations. However, 
the particular method used will vary from representation to representation. For exemplary purposes, the description 
below of the preferred embodiments, considers input images 109 as having been stored using the RGE components 
specified by the Society of Motion Picture and Television Engineers (SMPTE-C RGB). The white point is taken to be 
the standard CIE D65 white point, and the gamma is assumed to be either 1 (linear RGB) or 2.2 (gamma corrected 

20 for NTSC). The transformation and hence the contrast enhancement results will vary somewhat if different basis values 
are assumed, however it has been found that these assumptions work well for colour images that are designed to be 
viewed on a computer monitor. A more detailed discussion of these parameters is given in R.W.G. Hunt, The Repro- 
duction of Colour in Photography, Printing & Television, Fountain Press, Tolworth, England, 1987. 
[0055] Of course, a colour image 1 09 may have been stored in a format containing a luminance component, in which 

25 case the contrast enhancer 205 can skip the transformation step 501 . One embodiment uses the YCrCb colour space. 
Alternative embodiments, use CIELab, YUV, or YIQ colour spaces. 

[0056] In one embodiment, YCrCb is selected because it is a convenient representation for digital processing be- 
cause the chrominance components have easily computed maximum and minimum values. For example, if the RGB 
components each have an allowable range of 0 to 255, Y will fall between 0 and 255 as well, and the Cr and Cb values 
30 will both be between -1 27 and 1 28. A gamma corrected RGB image can be converted to YCrCb using a simple linear 
transformation, specifically, 



Y=.299R + .587G + .114B 

35 

Cr = .713 (R-Y) 



40 Cb = .564 (B-Y). 

[0057] The CIELab colour space, which was designed to be approximately perceptually uniform, is often used for 
digital image compression and transmission. The conversion to CIELab is somewhat more complicated, since the 
CIELab colour space is defined in terms of the standard CIE XYZ colour space. For gamma corrected RGB, this 
45 transformation is computed in three stages. First, convert to linear RGB. 

o = 255f — l 2 ' 2 G< = 255( — "f' 2 B, = 25$( — 

1 1.255 J 1 1 255 J 1 1 255 J 

50 

[0058] Next, a linear transformation is used to convert to CIE XYZ. 



X 




.3935 


.3653 


.1916 


y 




.2124 


.7011 


.0866 


z 




.0187 


.1119 


.9582 
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[0059] Finally, the data is converted to CIELab via the following equations. 



1 = 116 '(llo)- 16 

15 [ [ 100 ) { 108.89 J J 

where 

20 

f ( t) = ( C 1 t > ,008856 

7.7867t + .13793 t £ ,008856 

25 [0060] In one embodiment, these transformations are scaled and translated to produce integer values between 0 
and 255 for all colour components. In alternative embodiments, images may be digitised to other "bits-per-pixel" for- 
mats, e.g., four or sixteen. In those embodiments, the transformations are scaled to ranges appropriate for those 
formats. 

30 E. Luminance Histogram Clustering 

[0061] After the colour transformation, the data will consist of a luminance component 503 and two chrominance 
components 505 and 507. Next, the contrast enhancer 205 calls upon the histogram generator 209 to construct the 
luminance histogram 211 from the luminance data503 The histogram 211 is then input intoclusterizer213 to be divided 

35 into clusters, step 511. Assuming that luminance is recorded using one byte per pixel, the histogram 211 is a one 
dimensional array with 256 elements, so the required computation for clustering is independent of image size. 
[0062] The remapping of luminance values is based partially on cluster extent. Therefore, the exact threshold loca- 
tions can be quite important. The clusterizer 213 strives to divide the histogram such that flat regions are separate 
clusters, as illustrated in Figure 6. Figure 6 is an exemplary luminance histogram 600 having two peaks, one 601 in 

40 the low-luminance range and one 603 in the high-luminance range. These two peaks 601 and 603 are separated by 
a flat region 605 with low pixel counts. The clusterizer 213 strives to divide the exemplary histogram 600 such that flat- 
region 605 is a separate cluster. 

[0063] The clusterizer 213 uses a maximum likelihood technique for splitting clusters in conjunction with a tree grow- 
ing algorithm to determine the next cluster to split. Initially the clusterizer 213 views the histogram 211 as a single 
45 cluster, and successively splits existing clusters to generate the final result of a histogram 211 divided into several 
clusters. 

[0064] The procedure used by the clusterizer 213 is illustrated in Figure 7 and a flow-diagram of the procedure is 
shown in Figure 8. Initially the histogram 211 is viewed as one cluster 701, step 801. The cluster 701 is represented 
in a splitting tree 703 as node 701'. Next, the cluster 701 is split into two clusters 702a and 702b, which in turn are 

50 represented in the splitting tree 703 as nodes 702a' and 702b'. After several iterations, the histogram has been split 
into clusters 705a-e, corresponding to leaf-nodes 705a' through 705e' in the splitting tree 703. 
[0065] To decide which cluster to split, the clusterizer 213 searches among the existing clusters for the cluster that 
looks the least like a peak or a flat-region, step 803. To make that determination, the clusterizer 213 computes two 
different log likelihoods for each existing cluster, first assuming a Gaussian distribution and then a uniform distribution. 

55 a Gaussian distribution is taken to approximate a peak, and a uniform distribution is taken to approximate a flat-region. 
Therefore, a cluster that deviates a great deal from both a Gaussian distribution and a uniform distribution is a good 
candidate to split into additional clusters. 

[0066] The data values in each cluster are assumed to be independent and identically distributed (iid), and the pa- 
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rameters of the Gaussian distribution are taken to be the values that maximise the log likelihood. The log likelihood 
equations for a given cluster C under the Gaussian and uniform assumptions can thus be computed via (1) and (2), 
respectively. 

LL a = -i? (l+log(2BF 2 )) (1) 



LLy --/Vlog(|C|) (2) 

where 



* = fX 1 histii] 

Nice 



P 2 = <i-^) 2 histlil 

Nice 

\c\ =Ei 

i€C 

and the array value hist[i] is the number of image pixels in histogram bin i. 

[0067] The log likelihood associated with each cluster is then taken to be the maximum of (1) and (2). The cluster 
with minimum log likelihood is chosen as the next cluster to split. Thus, the algorithm chooses the existing cluster that 
looks the least like either a Gaussian or a uniform density. The average (per pixel) log likelihood of the chosen cluster 
is stored for each split. When this value is large enough in comparison with the value for the previous split, the clustering 
procedure terminates, step 804. The exact stopping criterion is given in Table 1 : pseudo-code for algorithm to divide 
histogram into clusters, 
initialize old11 = y, New1 1 = 0, and M = 1 
while (M < MaxM) 
for i = 1 to M 

use (1) to compute LL G (i) 
use (2) to compute LLy (i) 
set LL(i) to max(LL G (i), LL y (i)) 
set N(i) = number of pixels in cluster i 
setc = arg min. M M (LL(i)) 
setNew11 = LL(c)/N(c) 
if (New11 > 0.8*Old11) Stop 
reindex clusters C+1...M to C+2...M+1 

use (5) to split cluster c into 2 new clusters indexed c and c + 1 
if (c > 1) recompute boundary between c - 1 and c using (5) 
if (c < M) recompute boundary between c + 1 and c + 2 using (5) 
set M = M + 1 
setOld11 = New11 
terminate with M clusters. 
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Table 1 : Pseudo-code for algorithm to divide histogram into clusters. 

[0068] After selecting a cluster to split, the clusterizer 213 determines the point, within the range represented by the 
cluster, to split that cluster, step 805. To select the boundary point, the clusterizer 21 3 looks for uniformity in the histo- 
gram patterns on either side of the boundary point. For example, in the cluster 701 of Figure 7, the pixel counts in 
histogram bins left of the dividing line 707 are all relatively low and the pixel counts in the histogram bins to the right 
of dividing line 707 are relatively high. Therefore, the index at dividing line 707 is selected as the splitting boundary 
[0069] As a first step to select the boundary between the two new clusters, the clusterizer 213 models the data within 
the selected cluster as iid random samples from a mixture distribution. The distributions in the mixture are both discrete 
uniform distributions U (.,.), with the first extending from histogram bin 1 through bin m-l and the second from in through 
r-1, where the cluster extends from 1 through r-1 and in is between 1 and r. Each data value in the cluster is U(1,m) 
with some probability D and U(m,r) with probability (1 - D). The likelihood of the histogram data y in the cluster is thus 

GST- 

where 

N ± = ^hist[i] . 

[0070] To select the boundary at which to split the selected cluster, the clusterizer 213 seeks to pick m and d values 
such as to maximise (3). The value m will then be the boundary between the two new clusters. Taking the log of (3) gives 

losrp <j4m,D> = * t 1°9(-^i) + < n - n i )1os (t^) ' (4) 

[0071] The maximum likelihood estimate of D is computed by setting the partial derivative of (4) with respect to D to 
zero. 

This gives D = N 1 / N. Substituting into (4), selection of the boundary point is done by maximising the following function 
of m: 

logp<y|n,,D) - Wl log( 7? ^ IT ).(N-^ log ( J™^ ) (5 
= -fflog (r-I>+ ND |^|^Zlj , 

where 

D(p\\q) = p\og (g) + (1 - p) log (^j?) 

is the Kullback-Leibler distance between two binomial distributions with probabilities p and q, respectively (S. Kuliback 
and R.A. Leibler, On information and sufficiency. Annals of Mathematical Statistics, 22:79-86, 1951). Note that N 1 is a 
function of m. 

[0072] The clusterizer 213 maximises (5) using an exhaustive search on 1 +1 r-1 . This search gives the boundary 
between the two new clusters. 

[0073] Finally, the cluster boundaries to either side of the original cluster are recomputed given the result of the split, 
step 807. Table 1 is a pseudo-code listing of the clustering algorithm. 

[0074] In addition to terminating as described above in conjunction with Step 804, the procedure also terminates if 
the original histogram has been split into a maximum number of clusters, step 809. 
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F. Luminance Remapping 

[0075] After the clusterizer 213 has divided the luminance histogram into clusters, the contrast enhancer 205 calls 
upon the histogram luminance remapper 217 to remap the luminance values to new values based on cluster width, 
the number of pixels in each cluster, and the original histogram extent, step 513 of Figure 5. 
[0076] Let 



N(i) = number of pixels in cluster i 

Nt = total number of pixels in the image 

1 (i), r(i) = left and right limits of clusters i before modification 

a, b = left and right limits of histogram before modification 

a', b' = desired left and right limits of histogram after modification 



[0077] where a and b are defined to be the tightest bounds such that less than 1% of the image pixels fall to either 
side of the interval [a,...,b-1]. Given the desired values for a' and b', the luminance histogram remapper 217 remaps 
the cluster boundaries according to 



r . (i) =J . U) .j„(^).a-.,(^p)}<^', . ,7, 

where the parameter w is a weighting factor between 0 and 1 . The maximum dynamic range is achieved by setting a' 
and b' to the allowable limits, usually 0 and 256, respectively. However, most computer monitors and printers contain 
an offset below which all values map to black. The offset can vary from about 30 to 70 depending on the output device. 
In such applications, these parameters are set as follows: a' = 30, b' = 256 as the desired range. 

G. Histogram equalisation on clusters 

[0078] After the histogram remapper 217, using (6) and (7), has computed the new cluster boundaries, histogram 
equalisation is performed, step 51 5, by histogram equaliser/stretcher 21 9, on each cluster separately so the remapped 
data within each cluster is distributed as close to uniformly as possible, step 51 7. Alternatively, the histogram equaliser/ 
stretcher 219 can be used to stretch the data, so that the remapped data extends to the cluster boundaries but the 
shape of the histogram within a cluster remains unchanged. 

[0079] Figure 9 shows the results of performing this modification on a given image histogram. In this example, the 
remapping was performed using weighting factors of 0, 0.5, and 1 . The spikes that appear in the modified histograms 
are due to quantization effects during the pixel remapping, step 517. The data in any histogram bin can be remapped 
to a different value, but it cannot be split among multiple bins. For example, suppose bins 1 00 through 1 02 each contain 
1000 pixels initially. If we wish to remap this range to bins 110 through 113, we have 3 input bins being mapped to 4 
output bins. One of the output bins (1 1 0, 1 1 1 , 1 1 2, or 1 1 3) must therefore be empty while the other three contain 1 000 
pixels apiece. This will result in a downward spike in the output histogram. Similarly, if a given input range is mapped 
to a smaller output range, upward spikes will result. Both effects can be seen in Figure 9b-d. 

[0080] Note that a weighting factor of w = 0 in (7) will result in the relative cluster widths remaining unchanged after 
the remapping. This can limit the amount of contrast enhancement to too modest a level for images with clusters at 
both bright and dark luminance levels, as illustrated in Figure 9b. 

[0081] A weighting factor of w = 1, on the other hand, gives exactly the same mapping as histogram equalisation 
performed on the entire histogram; the clusters do not affect the mapping. As noted earlier, this can cause excessive 
luminance shifts for some images. 

[0082] In the preferred embodiment, a weighting factor of 0.5 is used, because this weighting factor has been found 
to result in good image quality over a large range of images. 
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H. Using Chrominance values to Prevent Overflow 

[0083] Thus far the remapping of pixel values has depended only on the luminance component. However, the max- 
imum allowable luminance value for a pixel depends on the chrominance of that pixel. If the luminance value exceeds 
this maximum, at least one of the colour components will be greater than the saturation level (usually 255) when the 
value is transformed back to RGB. To prevent this overflow from occurring, the contrast enhancer 205 calls upon a 
data checker 221 to use the chrominance values for each pixel to set a luminance threshold, step 51 9. If the remapped 
luminance exceeds this threshold for a given pixel, the data checker 221 sets the pixel luminance to the threshold value. 
[0084] The computation of a luminance threshold given the chrominance values is straightforward for the YCrCb 
colour space since Cr is proportional to R-Y and Cb is proportional to B-Y. For the YCrCb colour space, we can easily 
ensure that the red and blue components are not driven into saturation. In practice, colour shifts are very visible in skin 
tones, which are predominantly red, and in sky, for which blue dominates. However, saturation of the green component 
is generally less visible than for red and blue, so in an embodiment with reduced computations, green component 
overflow check is not necessary in the YCrCb implementation. 

[0085] For the CIELab colour space the overflow checking is less straightforward. There is no simple way to relate 
the overflow of a specific RGB component to Lab values. Therefore, in a CIELab embodiment, a two dimensional 
lookup table, indexed by the quantized a and b values, that contains the luminance threshold, is employed. To construct 
the look-up table, the eight corners of RGB space (RGB = (0,0,0), (0,0,255), (0,255,0), etc.) were first converted to 
CIELab colour space. The corner values constitute eight values for the lookup table. The table is filled in using bilinear 
interpolation between these corner values. The resulting thresholds are not precisely correct, but any overflow errors 
that still occur are small enough that only minor colour shifts result. Thus, visual degradation is kept small. This pro- 
cedure requires a 256 x 256 lookup table of unsigned characters, but its operation is very fast. Alternatively, the bilinear 
interpolation could be performed for each image pixel in turn, which would not require the extra memory but would be 
much slower. 

I. Convert to RGB colour space 

[0086] Finally, having computed remapped values for all pixels in the image 109, the image is converted back to 
RGB space, step 521. The resulting image may then be stored back into RAM 105 or written onto a hard disk 119. 
Alternatively, the image may be printed on a printer 115 or output on another output device. 

J. Conclusion 

[0087] The preferred contrast enhancer 205 may be used to improve the image quality of an image at any stage in 
the process flow of Figure 3. For example, an input device, such as a digital camera may include a contrast enhancer 
205 and improve the contrast of an image prior to transferring the image to a computer 1 01 . Alternatively, the computer 
1 01 may include the contrast enhancer 205 to improve image quality of images prior to printing or displaying the images. 
The contrast enhancer 205 may also be incorporated into image processing software packages such as computerised 
"dark-rooms". Yet another alternative use of a contrast enhancer 205 is as part of a printer driver 125. A further alter- 
native is to place the contrast enhancer 205 "on-board" an output device, e.g., a printer 115 in its internal control 
mechanism. 



Claims 

1. A method of operating a computer to produce contrast enhanced images, wherein an image is composed of a 
plurality of pixels, comprising the steps of: 

inputting (1 09) in digital form an original image into a storage unit of said computer; 

producing (509) a histogram of said digital form having a first axis corresponding to a first measurable property 
and asecond axis corresponding to acount of pixels having a particular value for said first measurable property; 
detecting patterns in the histogram by detecting regions in said histogram that resemble Gaussian distributions 
and detecting other regions in said histogram that resemble uniform distributions and dividing (511) the his- 
togram into clusters based on said patterns; 

performing (515) histogram equalisation on each cluster separately thereby producing a modified histogram; 
and 

using (517) said modified histogram to adjust the value of said first measurable property in said digital form, 
thereby producing said contrast enhanced image. 
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2. A method as in claim 1, wherein said dividing step comprises the step of partitioning the histogram into clusters 
of equal size. 

3. A method as in claim 1 or 2, wherein for an RGB image the method comprises the step of converting said RGB 
image into a format having a luminance component. 

4. A method as in any preceding claim, wherein said first measurable property is a light intensity measure and wherein 
for each light intensity value said count of pixels is a count of pixels having that particular light intensity, and wherein 
said histogram has an extent equal to the highest and lowest light intensity values in said original image, each 
cluster representing a range in said histogram, the sum of the histogram values for the range corresponding to 
the number of pixels having a light intensity within the range of the cluster; the method comprising the step of 
remapping cluster boundaries proportionally to the cluster's range (Ro), said number of pixels within said range 
(Ni), the total number of pixels in the image (Nt), and said histogram extent (Eo). 

5. A method as in claim 4, wherein said remapping step comprises the step of selecting a desired range for the 
histogram (Ed) of said output image; and stretching the width of each cluster according to the relationship: 













w 




+ a-w) 








i E o) 
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wherein w is a weighting factor ranging from 0 to 1 .0 and R n is the new cluster width. 

6. A method as in any preceding claim, wherein for an image with a digital form which comprises a second measurable 
property, said generating step comprises the step of for each pixel in said contrast enhanced image, limiting the 
value of said first measurable property as a function of the value of said second measurable property. 

7. A method as in claim 6, comprising the steps of generating a table of maximum values of said first measurable 
property as a function of said second measurable property; said limiting step comprising, for each value V of said 
second measurable property, the step of retrieving from said table the maximum value M for said first measurable 
property corresponding to said value V and if said value V exceeds said maximum value H setting the value V to H. 

8. An image processing system for producing contrast enhanced images, wherein an image is composed of a plurality 
of pixels, comprising: 

a storage unit (105) operable to store a digital form of an original image; 

a histogram producing unit (209) operable to produce a histogram of a said digital form having a first axis 
corresponding to a first measurable property and a second axis corresponding to a count of pixels having a 
particular value for said first measurable property; 

pattern detecting means operable to detect patterns in a said histogram by detecting regions in the histogram 
that resemble Gaussian distributions and detecting other regions in the histogram that resemble uniform dis- 
tributions and to divide (511) the histogram into clusters based on the detected patterns; 
a histogram equaliser (21 9) operable to perform histogram equalisation on each cluster separately thereby to 
produce a modified histogram; and 

adjustment means operable by using a said modified histogram to adjust the value of said first measurable 
property in said digital form, thereby to produce a contrast enhanced image. 

9. An image processing system as in claim 8, wherein said first measurable property is a light intensity measure and 
wherein for each light intensity value said count of pixels is a count of pixels having that particular light intensity, 
and wherein a said histogram has an extent equal to the highest and lowest light intensity values in the original 
image, each cluster representing a range in said histogram, the sum of the histogram values for the range corre- 
sponding to the number of pixels having a light intensity within the range of the cluster; the histogram equaliser 
comprising the means for remapping cluster boundaries proportionally to the cluster's range (ro), the number of 
pixels within the range (Ni), the total number of pixels in the image (Nt), and histogram extent (Eo). 

10. An image processing system as in claim 9, wherein said remapping means is operable to select a desired range 
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for the histogram (Ed) of the output image and to stretch the width of each cluster according to the relationship: 



wherein w is a weighting factor ranging from 0 to 1 .0 and R n is the new cluster width. 

An image processing system as in claim 8, 9 or 1 0, comprising an image translator (207) containing transformations 
between a first image format and a second image format that includes said first property. 

An image processing system as in any one of claims 8 to 11, comprising an image generator connected to said 
histogram data structure and said original image, and operable to generate a contrast enhanced image from said 
histogram data structure and said original image. 



Ein Verfahren zum Betreiben eines Computers, urn Kontrast-erhohte Bilder zu erzeugen, bei dem ein Bild aus 
einer Mehrzahl von Pixeln besteht, mit folgenden Schritten: 

Eingeben (109) in digitaler Form eines Originalbilds in eine Speichereinheit des Computers; 

Erzeugen (509) eines Histogramms der digitalen Form, das eine erste Achse, die einer ersten meGbaren 
Eigenschaft entspricht, und eine zweite Achse aufweist, die einem Zahlwert von Pixeln entspricht, die einen 
bestimmten Wert fur die erste meGbare Eigenschaft aufweisen; 

Erfassen von Mustern in dem Histogramm durch ein Erfassen von Bereichen in dem Histogramm, die GauG' 
schen Verteilungen gleichen, und Erfassen von anderen Bereichen in dem Histogramm, die gleichformigen 
Verteilungen gleichen, und Teilen (51 1 ) des Histogramms in Bundel auf der Basis der Muster; 

Durchfuhren (515) eines Histogrammabgleichs fur jedes Bundel getrennt, wodurch ein modifiziertes Histo- 
gramm erzeugt wird; und 

Verwenden (51 7) des modifizierten Histogramms, urn den Wert der ersten meGbaren Eigenschaft in der digi- 
talen Form einzustellen, wodurch das Kontrast-erhohte Bild erzeugt wird. 

Ein Verfahren nach Anspruch 1 , bei dem der Schritt eines Teilens den Schritt eines Partitionierensdes Histogramms 
in Bundel mit gleicher GroGe aufweist. 

Ein Verfahren nach Anspruch 1 oder 2, bei dem das Verfahren fur ein RGB-Bild den Schritt eines Umwandelns 
des RGB-Bilds in ein Format aufweist, das eine Luminanzkomponente aufweist. 

Ein Verfahren nach einem der vorhergehenden Anspruche, bei dem die erste meGbare Eigenschaft ein Lichtin- 
tensitatsmaG ist, und bei dem fur jeden Lichtintensitatswert der Zahlwert von Pixeln ein Zahlwert von Pixeln ist, 
die diese bestimmte Lichtintensitat aufweisen, und bei dem das Histogramm eine Ausdehnung aufweist, die gleich 
dem hochsten und niedrigsten Lichtintensitatswert in dem Originalbild ist, wobei jedes Bundel einen Bereich in 
dem Histogramm darstellt, wobei die Summe der Histogrammwerte fur den Bereich der Anzahl von Pixeln ent- 
spricht, die eine Lichtintensitat in dem Bereich des Bundels aufweisen; 

wobei das Verfahren den Schritt eines erneuten Abbildens von Bundelbegrenzungen proportional zu dem Bun- 
delbereich (Ro), der Anzahl von Pixeln in dem Bereich (Ni), der Gesamtanzahl von Pixeln in dem Bild (Nt), und 
der Histogrammausdehnung (Eo) aufweist. 

Ein Verfahren nach Anspruch 4, bei dem der Schritt eines erneuten Abbildens den Schritt eines Auswahlens eines 
gewunschten Bereichs fur das Histogramm (Ed) des Ausgabebilds und eines Stretchens der Breite jedes Bundels 
gemaG der Beziehung umfaGt: 



5 




13 
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wobei w ein Gewichtungsfaktor ist, der in einem Bereich von 0 bis 1,0 liegt, und R n die neue Bundelbreite ist. 

6. Ein Verfahren nach einem der vorhergehenden Anspruche, bei dem fur ein Bild mit einer digitalen Form, das eine 
zweite meGbare Eigenschaft aufweist, der Schritt eines Erzeugens den Schritt eines Begrenzens des Werts der 

10 ersten meGbaren Eigenschaft als eine Funktion des Werts der zweiten meGbaren Eigenschaft fur jedes Pixel in 

dem Kontrast-erhohten Bild aufweist. 

7. Ein Verfahren nach Anspruch 6, das die Schritte eines Erzeugens einer Tabelle von maximalen Werten der ersten 
meGbaren Eigenschaft als eine Funktion der zweiten meGbaren Eigenschaft aufweist; 

15 wobei der Schritt eines Begrenzens fur jeden Wert V der zweiten meGbaren Eigenschaft den Schritt eines Wie- 

dergewinnens des maximalen Werts M fur die erste meGbare Eigenschaft, der dem Wert V entspricht, aus der 
Tabelle und eines Einstellens des Werts V auf H, wenn der Wert V den maximalen Wert H uberschreitet, aufweist. 



8. Ein Bildverarbeitungssystem zum Erzeugen von Kontrast-erhohten Bildern, bei dem ein Bild aus einer Mehrzahl 
20 von Pixeln besteht, mit folgenden Merkmalen: 

einer Speichereinheit (1 05), die betrieben werden kann, urn eine digitale Form eines Originalbilds zu speichern; 

einer Einheit (209) zum Erzeugen eines Histogramms, die betrieben werden kann, urn ein Histogramm der 
25 digitalen Form zu erzeugen, das eine erste Achse, die einer ersten meGbaren Eigenschaft entspricht, und 

eine zweite Achse aufweist, die einem Zahlwert von Pixeln entspricht, die fur die erste meGbare Eigenschaft 
einen bestimmten Wert aufweisen; 



einer Mustererfassungseinrichtung, die betrieben werden kann, urn Muster in dem Histogramm durch ein Er- 
30 fassen von Bereichen in dem Histogramm, die GauG'schen Verteilungen gleichen, und ein Erfassen anderer 

Bereiche in dem Histogramm, die gleichformigen Verteilungen gleichen, zu erfassen und das Histogramm auf 
der Basis der erfaGten Muster in Bundel zu teilen (511); 

einem Histogrammabgleicher (21 9), der betrieben werden kann, urn einen Histogrammabgleich fur jedes Bun- 
35 del getrennt durchzufuhren, um dadurch ein modifiziertes Histogramm zu erzeugen; und 



einer Einstelleinrichtung, die unter Verwendung des modifizierten Histogramms betrieben werden kann, um 
den Wert der ersten meGbaren Eigenschaft in der digitalen Form einzustellen, um dadurch ein Kontrasterhoh- 
tes Bild zu erzeugen. 

40 

9. Ein Bildverarbeitungssystem nach Anspruch 8, bei dem die erste meGbare Eigenschaft ein LichtintensitatsmaG 
ist, und bei dem fur jeden Lichtintensitatswert der Zahlwert von Pixeln ein Zahlwert von Pixeln ist, die diese be- 
stimmte Lichtintensitat aufweisen, und bei dem das Histogramm eine Ausdehnung aufweist, die gleich dem hoch- 
sten und niedrigsten Lichtintensitatswert in dem Originalbild ist, wobei jedes Bundel einen Bereich in dem Histo- 

45 gramm darstellt, wobei die Summe der Histogrammwerte fur den Bereich der Anzahl von Pixeln entspricht, die 

eine Lichtintensitat in dem Bereich des Bundels aufweisen; 

wobei der Histogrammabgleicher die Einrichtung zum erneuten Abbilden von Bundelbegrenzungen proportional 
zu dem Bundelbereich (Ro), der Anzahl von Pixeln in dem Bereich (Ni), der Gesamtanzahl von Pixeln in dem Bild 
(Nt), und der Histogrammausdehnung (Eo) aufweist. 

50 

10. Ein Bildverarbeitungssystem nach Anspruch 9, bei dem die Einrichtung zum erneuten Abbilden betrieben werden 
kann, um einen gewunschten Bereich fur das Histogramm (Ed) des Ausgabebilds auszuwahlen und die Breite von 
jedem Bundel gemaG der Beziehung: 



14 
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zu stretchen, wobei w ein Gewichtungsfaktor ist, der in einem Bereich von 0 bis 1,0 liegt, und R n die neue Bun- 
delbreite ist. 

11. Ein Bildverarbeitungssystem nach einem der Anspruche 8 bis 10, das einen Bildubersetzer (207) aufweist, der 
Transformationen zwischen einem ersten Bildformat und einem zweiten Bildformat, das die erste Eigenschaft 
umfaBt, enthalt. 

12. Ein Bildverarbeitungssystem nach einem der Anspruche 8 bis 11 , das einen Bilderzeuger, der mitder Histogramm- 
datenstruktur und dem Originalbild verbunden ist, und der betrieben werden kann, urn ein Kontrasterhohtes Bild 
aus der Histogrammdatenstruktur und dem Originalbild zu erzeugen. 



Revendications 

1 . Un procede de mise en oeuvre d'un ordinateur pour reproduire des images a contraste rehausse, dans lequel une 
image est composee d'une serie de pixels, qui comprend les etapes consistant a: 

entrer (109) sous forme numerique une image d'origine dans une unite de memoire dudit ordinateur; 
produire (509) un histogramme de ladite forme numerique qui comporte un premier axe correspondant a une 
premiere propriete mesurable et un deuxieme axe correspondant a un compte de pixels d'une valeur particu- 
liere pour ladite propriete mesurable; 

detecter des configurations dans I'histogramme en detectant des regions dudit histogramme qui ressemblent 
a des repartitions gaussiennes et en detectant d'autres regions dudit histogramme qui ressemblent a des 
repartitions uniformes et diviser (511 ) I'histogramme en groupes sur la base desdites configurations; 
effectuer (515) une egalisation d'histogramme separement sur chaque groupe, en produisant ainsi un histo- 
gramme modifie; et 

utiliser (517) ledit histogramme modifie pour ajuster la valeur de ladite premiere propriete mesurable dans 
ladite forme numerique, en produisant ainsi une image a contraste rehausse. 

2. Un procede selon la revendication 1 , dans lequel ladite etape de division comprend I'etape consistant a segmenter 
rhistogramme en groupes d'egale dimension. 

3. Un procede selon la revendication 1 ou 2, dans lequel le procede comprend, pour une image RGB, I'etape con- 
sistant a convertir ladite image RGB en un format a composante de luminance. 

4. Un procede selon I'une quelconque des revendications precedentes, dans lequel ladite premiere propriete mesu- 
rable est une mesure d'intensite lumineuse et dans lequel ledit compte de pixels pour chaque valeur d'intensite 
lumineuse est un compte des pixels de cette intensite lumineuse particuliere, et dans lequel I'etendue de I'histo- 
gramme est egale aux valeurs maximale et minimale d'intensite lumineuse dans ladite image d'origine, chaque 
groupe representant une plage dudit histogramme, la somme des valeurs de I'histogramme pour la plage corres- 
pondant au nombre de pixels a intensite lumineuse comprise a I'interieur de la plage du groupe; le procede com- 
prenant I'etape consistant a reconfigurer des limites de groupe proportionnellement a la plage du groupe (Ro), 
audit nombre de pixels a I'interieur de ladite plage (Ni), au nombre total de pixels de I'image (Nt) et a ladite etendue 
(Eo) d'histogramme. 

5. Un procede selon la revendication 4, dans lequel ladite etape de reconfiguration comprend I'etape consistant a 
selectionner une plage souhaitee pour I'histogramme (Ed) de ladite image de sortie; et a etaler la largeur de chaque 
groupe selon la relation: 
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ou w est un facteur de ponderation compris entre 0 et 1 ,0 et R n est la nouvelle largeur de groupe. 



Un procede selon Tune quelconque des revendications precedentes, dans lequel ladite etape de generation com- 
prend, pour une image a forme numerique qui comprend une deuxieme propriete mesurable, I'etape consistant a 
limiter pour chaque pixel de ladite image a contraste rehausse la valeur de ladite premiere propriete mesurable 
en fonction de la valeur de ladite deuxieme propriete mesurable. 



7. Un procede selon la revendication 6, qui comprend les etapes consistant a engendrer une table de valeurs maxi- 
males de ladite premiere propriete mesurable en fonction de ladite deuxieme propriete mesurable; ladite etape 
de limitation comprenant, pour chaque valeur V de ladite deuxieme propriete mesurable, I'etape consistant a res- 
tituer de ladite table la valeur maximale M de ladite premiere valeur mesurable correspondant a ladite valeur V et 
a regler a H la valeur V si ladite valeur V depasse ladite valeur maximale H. 

8. Un systeme de traitement d'images pour produire des images a contraste rehausse, dans lequel une image est 
composee d'une serie de pixels, comprenant: 

une unite de memorisation (105) qui peut intervenir pour memoriser une forme numerique d'une image d'ori- 
gine; 

une unite de production (209) d'histogrammes qui peut intervenir pour produire un histogramme d'une dite 
forme numerique qui comporte un premier axe correspondant a une premiere propriete mesurable et un 
deuxieme axe correspondant a un compte de pixels d'une valeur particuliere pour ladite propriete mesurable; 
un moyen de detecter de configurations qui peut intervenir pour detecter des configurations dans I'histogram- 
me en detectant des regions de I'histogramme qui ressemblent a des repartitions gaussiennes et en detectant 
d'autres regions de I'histogramme qui ressemblent a des repartitions uniformes et diviser (511 ) I'histogramme 
en groupes sur la base desdites configurations; 

un egalisateur d'histogramme (21 9) qui peut intervenir pour effectuer une egalisation d'histogramme separe- 
ment sur chaque groupe, afin de produire ainsi un histogramme modifie; et 

un moyen d'ajustement susceptible d'intervenir en utilisant un dit histogramme modifie pour ajuster la valeur 
de ladite premiere propriete mesurable dans ladite forme numerique, afin de produire ainsi une image a con- 
traste rehausse. 



9. Un systeme de traitement d'images selon la revendication 8, dans lequel ladite premiere propriete mesurable est 
une mesure d'intensite lumineuse et dans lequel ledit compte de pixels pour chaque valeur d'intensite lumineuse 
est un compte des pixels de cette intensite lumineuse particuliere, et dans lequel I'etendue d'un dit histogramme 
est egale aux valeurs maximale et minimale d'intensite lumineuse dans I'image d'origine, chaque groupe repre- 
sentant une plage dudit histogramme, la somme des valeurs de I'histogramme pour la plage correspondant au 
nombre de pixels a intensite lumineuse comprise a I'interieur de la plage du groupe; I'egalisateur d'histogramme 
comprenant le moyen de reconfiguration de limites de groupe proportionnellement a la plage (Ro) du groupe, au 
nombre de pixels a I'interieur de la plage (Ni), au nombre total de pixels de I'image (Nt) et a I' etendue (Eo) d'his- 
togramme. 



10. Un systeme de traitement d'images selon la revendication 9, dans lequel ledit moyen de reconfiguration peut 
intervenir pour selectionner une plage souhaitee pour 1'histogramme (Ed) de I'image de sortie et pour etaler la 
largeur de chaque groupe selon la relation: 



Rn = 
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ou w est un facteur de ponderation compris entre 0 et 1 ,0 et R n est la nouvelle largeur de groupe. 
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Un systeme de traitement d'images selon la revendication 8, 9 ou 10, comprenant un traducteur (207) d'images 
qui contient des transformations entre un premier format d'image et un deuxieme format d'image qui inclut ladite 
premiere propriete. 

Un systeme de traitement d'images selon Tune quelconque des revendications 8 a 11 , qui comprend un generateur 
d'images connecte a ladite structure de donnees d'histogrammes et a ladite image d'origine, et peut intervenir 
pour engendrer une image a contraste rehausse a partir de ladite structure de donnees d'histogramme et de ladite 
image d'origine. 
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